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[57] ABSTRACT 

A method is detailed to determine the density of parti- 
cles wherein the closing velocity is known between the 
impacting particles and a plate of known material. Ei- 
ther the shockwave velocity or the material velocity 
produced in the plate upon impact by an unknown 
material particle is determined and compared with the 
corresponding shockwave or material velocity that 
would be produced by different known material parti- 
cles having the same closing velocity upon impact with 
the plate. The unknown material particle density is de- 
rived by obtaining a coincidence of the shockwave ve- 
locity or material velocity conditions initially, pro- 
duced upon impact between the known material plate 
and one of the different material particles and from 
the fact that shockwave velocity and material velocity 
are ordered on the impacting particle material density 
alone. 

2 Claims, 19 Drawing Figures 
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DETERMINING PARTICLE DENSITY USING 
KNOWN MATERIAL HUGENIOT CURVES 

ORIGIN OF THE INVENTION 

The invention described herein was made by an em- 
ployee of the United States Government and may be 
manufactured and used by or for the Government for 
governmental purposes without the payment of any 
royalties thereon or therefor. 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The invention relates to the determination of the 
density of particulate matter, and has particular use 
with spacecraft. 

2. Description of the Prior Art 

It is often necessary to determine the constituents 
and density of the individual particles, in space, in addi- 
tion to the mass distribution of particulate matter. 
Properly designed spacecraft meteoroid shielding uti- 
lizing the meteoroid bumper concept depends inter alia 
on the density of meteoroid particles. 

Any practicable device used to determine the density 
of particles must be simple, lightweight and compact in 
order to be placed on small unmanned exploratory 
spacecraft. Prior art devices are based on either captur- 
ing the particle and conducting an examination of it, or 
using a spectrograph to determine the particle constitu- 
ents from an impact flash generated when a particle 
collides with the spacecraft. 

SUMMARY OF THE INVENTION 

The instant application relates to the determination 
of the density of all known velocity cosmological mate- 
rial. There is at this time no prior art, previous methods 
or devices performing the function of the instant appli- 
cation. 

The invention provides a simple, compact and light- 
weight device because only one velocity measurement 
must be made. This velocity measurement can be made 
by recording impact flashes which impinge on a photo- 
diode. In addition, the instant invention can be fully in- 
tegrated with any capture device, and is also com- 
pletely complementary to any spectrograph system as 
outlined below, in the following sense. With a spectro- 
graph the elements in the particle and their relative 
amount are determined. Calculation, assuming a solid 
homogeneous particle, would then yield the particle 
density. In the instant invention, density is determined 
and calculation used to determine its principal constit- 
uents. For particles containing voids produced, for ex- 
ample, in the accreted particles in the asteroid belt, si- 
multaneous measurements of the particle constituents 
by spectrograph, and the particles density by devices 
described here, are needed to have an adequate de- 
scription of the impacting particle. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows the impending impact between a parti- 
cle having a known velocity and a plate of finite thick- 
ness; 

FIG. 2 shows the conditions present in the particle 
and plate of FIG. 1, shortly after impact; 

FIG. 3 illustrates a graphical procedure for determin- 
ing the pressure and material velocity of particle and- 
plate during impact; 
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FIG. 4 shows conditions reached in a known material 
plate impacted by an unknown material particle having 
a known closing velocity; 

FIG. 5 shows the use of trial Hugoniots in a shifted 
5 and reflected technique to determine an unknown par- 
ticle density; 

FIG. 6 are Hugoniot curves for a plurality of materi- 
als; 

FIGS, la and lb illustrate a means to measure plate 
10 shockwave velocity; 

FIG. 8 shows the conditions of a finite thickness 
plate, with a reflected expansion wave from plate rear 
surface; 

FIG. 9 shows the impact pressure and particle veloc- 
15 ity produced by high speed impact and maximum spray 
velocity produced for the case > pta, 

FIGS. lOu and 10f> show a means to measure debris 
cloud velocity; 

FIG. 11 is a graph of light level versus time for the 
20 photodiode shown in RGS. 10a and lOi, when p^o > 
Plot 

FIG. 12 shows impact pressure, particle velocity and 
maximum spray density produced for hypervelocity im- 
pact when ppo = Pto'y 

25 fig. 13 is a graph of light level versus time for the 
photodiode shown in FIGS. lOn and lOh, when Pp„ = 
Ptor 

FIG. 14 shows impact pressures and plate velocities 
produced during successive shockwave and expansion 
30 vvave passages, when Pp„ < pm\ 

FIG. 15 is a graph of light level versus time for the 
photodiode shown in FIGS. lOn and 101), when fpo < 
Pto, 

FIG. 16 is a graph showing the intersection of the be- 
33 ryllium Hugoniot with the Hugoniots of various other 
materials; and 

FIG. 17 is a table of the material and free surface ve- 
locities produced in a beryllium plate impacted by par- 
ticles of various materials. 

40 

DETAILED DESCRIPTION OF THE INVENTION 

In the method of the invention, density is determined 
and calculations are made to determine the principal 
constituents. One restriction is that the particle closing 
velocity must be less than 2 cm/microsec., in order that 
particle and spacecraft surface not vaporize when one 
of the devices proposed here is utilized. In addition, 
one other restriction is that the closing velocity be 
known. The closing velocity is a known quantity in cer- 
tain applications, as for example, in determining parti- 
cle density in the heads of comets such as the d’Arrest 
comet, for particles in the asteroid belt, and particles 
in orbit about such planets as Jupiter and Saturn. If the 
particle-spacecraft closing velocity is known, the angle 
^ of particle impact with a spacecraft surface is also 
known and can be controlled. 

SYMBOLS 

C — isentropic speed of sound 
W — void width between a thin plate and a backstop 

F — impact flash 

_ J — right circular cylinder length 
L — characteristic length of meteoroid 
P — pressure 
R — expansion wave 
S — shockwave 



3 , 842,656 


h — plate thickness 
T — time 

u — material velocity 

PiQj Cpy Bpt Ppot ^ pi ( [OtPto “h ^pPpo ”1~ 

2BePpok'p] + 

((a,p,„ + flpPpo + 2BppppVp)'^ + 4(B,p,„ - BpPpp) 

(OpPppVp + Bp^poKpOl''') (2(B,p,p - B^pp)] ; (Bip,„ - 
BpAo) »* 0 

U — shockwave velocity 

Vp — particle velocity normal to plate surface 

P — density 

w — right circular cylinder width 
W — characteristic width of meteoroid 
X — position along impact axis 
a — material constant (cm/microsec.) 

B — material constant (dimensionless) 

E — characteristic meteoroid surface roughness 
length 

AT — time after impact 
6 — angle 

Subscripts 

0 — state of material originally 

1 — state of material after first shockwave passage or 
first shockwave and expansion wave passage 

2 — state of material after second shockwave passage 
or second shockwave and expansion wave passage 

3 — state of material after third shockwave passage 
or third shockwave and expansion wave passage 

p — particle 
t — plate 

tr — pertaining to expansion wave from plate rear 
surface 

fs — pertaining to plate rear surface 
max — pertaining to a maximum 
min — pertaining to a minimum 
It is necessary to consider in detail the actual colli- 
sion process which occurs when a hypervelocity parti- 
cle impacts a finite thickness target. Because the angle 
at which a particle impacts a spacecraft surface can be 
control led, o nly nora^, or perpendicular, impacts will 
be considered as occurring. FIG. 1 shows the condi- 
tions present prior to impact of a known velocity parti- 
cle 10 and a finite thickness plate 11. 

The particle has characteristic length L, width W, 
surface roughness length E, and the plate has thickness 
h. Consideration will be limited to impact cases such 
that , 

E < h < min (L,W). 

Upon particle-thin plate impact, a strong shockwave, 
Si, is driven into the plate as shown in FIG. 2 and an- 
other shockwave Sp is driven back into the particle. 
The momentum equations (2nd Rankine-Hugoniot re- 
lation) across the plate and particle shockwaves written 
relative to the undisturbed plate and particle materials 
are 

F(/ ~ PtoUtUn 

Ppl PptfUpMp/ 

( 2 ) 


to 


20 


30 


35 


40 


45 


where o subscripts denote the undisturbed ma terial 
state before shock passage and 1 subscripts deriotelhe 
material state after shock passage. In the shaded area 
of FIG. 2, P,/ and Pp/ are the . initially induced impact 
pressures in the plate and particle material, p,o and Pp„ 
are the original plate and particle material densities, 
and u,i with and Up, with U„ are the initially induced 
compressed material and shockwave velocity for the 
plate and particle relative to the undisturbed plate and 
particle material. For solids there exists a linear rela- 
tion between the shock velocity U and compressed ma- 
terial velocity u of the form 


U=a+Bu 


15 


( 3 ) 


where a(cm/microsec.) and B (dimensionless) are em- 
pirical constants for each material. Substitution of 
equation (3) into equations (1) and (2) yields two 
equations 


25 


P,i=Pio(a, -t- B,u„) u,t 


Ppl = Pi»(flp + BpUpi) Upt. 


( 4 ) 


( 5 ) 


To connect the pressure and material velocity in each 
material in the shaded region of FIG. 2, it may be noted 
that across any part of the particle-plate interface, for 
example, A-B as we are interested in the initially in- 
duced impact states, Newton’s Second Law holds and 
there is no interpenetration of one material into an- 
other. This implies that 


and 


Pti — Ppl 


u,i = Vp- Upi 


(6) 


( 7 ) 


Equating equations (4) and (5) through equation (6) 

and noting from equation (7) 

that 


Upt =Vp- u,i 


50 


(8) 


55 


60 


one may solve for as a function of impact velocity 
if both particle and plate materials are known as fol- 
lows; 

pioia, + B,u„)u„ = ppp(ap + Bp[ Vp - u„] ) (Vp- u„) 
which leads to the quadratic 

(Bp,o - B^pp) 

"F (uiPto "F QpPpp "F 2BpPppVp) Uft 
(UpPppVp -F BpPppVp^) = 0 

which may be solved taking the plus sign of the qua- 
dratic as 


Utl = 


— latPto+ayPpo+2BpPp„Vp]+l(atPto+apPvo+2BpPppVp)‘+4(BtPto—BpPpp)(apPppVp+BpPppVp^)]^ 


2(BiPip BpPpo} 


= Vp—Upi; (B,pip-BpPpp)p^0 
(9) 
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Using equations (9) and (4), P,/ as a function of Vp 
may also be'solved for with the result 

Ptl Ppl Pio 4" SlUif (fl|, ptoy ^lly Pliy Ppo, 
u„ (a„ B„ p,p, Op, Bp, ppo, Pp). 


( 10 ) 


A straightforward method of determining U, would 
be to record the time, T, it takes the plate shockwave, 
S/, to propagate through a certain thickness plate h. 
TTien 


u, = hrr. 


( 11 ) 


With equations (9) and ( 10) one can calculate the ini- 
tial pressure and compressed material velocity pro- 
duced, knowing the particle and plate materials and the 
closing velocity. 

There is, however, a graphical procedure which may 
be used to determine the pressure and compressed ma- 
terial velocity produced. First, Hugoniot curves [equa- 
tions of the same form as equation (4)] are calculated 
for the plate and particle materials. On a plot (FIG. 3) 
which has pressure on the ordinate and material veloc- 
ity on the abscissa, the Hugoniot curve P„ = P,i(u,i) is 
plotted starting at «« = 0, P = 0. Starting at the impact 
velocity Vp on the abscissa and at P = 0 , the particle 
Hugoniot is plotted and reflected back through the 
plate Hugoniot. The intersection of the two Hugoniots 
gives the pressure and material velocity initially pro- 
duced in the particle and plate through impact. Thus 
the plot says an element of plate material initially at 
rest and at zero pressure is brought through action of 
a shockwave to some P,/ and u,i on the plate Hugoniot. 
Also, an element of particle material at initial velocity 
Vp and zero pressure is brought through action of a 
shockwave to some Pp/ and Fp - «p,. Equations ( 6 ) and 
(7), which specify that P</ and Pp/ must be the same and 
u,i= Fp - Up,, are precisely the conditions present at the 
intersection of the normal plate material Hugoniot and 
the shifted and reflected particle material Hugoniot. 

Now consider the situation as depicted in FIG. 4, as- 
suming the pressure P„ and compressed material veloc- 
ity u,i which has been produced in a known material 
plate by impact is known. Assume further that the par- 
ticle impact velocity is known but the particle material 
is unknown. To determine the unknown particle den- 
sity, trial known material Hugoniots shifted to Vp are 
plotted as shown in FIG. 5. 

The trial Hugoniot which passes through the known 
P„ and u,i on the plate Hugoniot will give the density of 
the unknown particle material because it is well known 
that the Hugoniots are ordered with density alone. FIG. 
6 , which was obtained from extrapolation of low pres- 
sure data supplied in “Shock Waves, High Pressure 
Physics and Chemistry,” Academic Press, pp. 209-291 
(1963), by Duvall and Fowles, illustrates this fact. The 
Hugoniot curves from gold through water may be seen 
to be ordered on density alone, and various room tem- 
perature properties of the materials such as ductility, 
yield strength, etc., have no effect. 

The problem now is to outline how the pressure P« 
and compressed material velocity u„ are to be deter- 
mined. Examination of equation ( 1 ) reveals that two 
quantities, U„ the shockwave velocity in the plate, and 
u,i, the compressed material velocity in the plate, are 
available for measurement. Equation (3) then reveals, 
that a measurement of either allows the other, as well 
us P„ to be calculated. 


FIG. 7 shows a schematic of a possible device to carry 
10 out this measurement. 

A photodiode structure 20 is used to record the im- 
pact flash 21 on the target front surface 22 and to start 
a time counter 23. By backing the plate with argon gas 
24, a strong shockwave 25 would be driven into the gas 
when the target shockwave reaches the target rear sur- 
face. Strongly shocked argon gas 26 glows brightly and 
this flash can be sensed by photodiode structure 20 to 
turn the time counter 23 off. Although the technique 
is relatively straight forward, it is severely restricted by 
the expansion waves 27 which originate at the periph- 
ery of the particle-plate contact surface. Starting at the 
periphery of the particle-plate contact surface, these 
expansion waves decompress the compressed particle 
and plate materials, and overtake and reduce the initial 
shockwave velocity in the particle and plate. Any meas- 
urement of shockwave velocity must be accomplished 
in a plate whose thickness is less than the depth at 
which the total shockfront suffers attenuation. In gen- 
30 eral the counter 23 must operate on the 1 X 1 to 1 
X 10"'“ sec. time scale. If this time measurement can be 
made, there should be no restriction on the density of 
meteoroids which can be measured with perhaps the 
sole exceptions of those characterized by a snowflake- 
35 type structure. 

Measurement of u„ 


To determine m„, the shockwave S, must again reach 
the plate rear surface unattenuated by expansion waves 
from the side and rear of the particle. Again, the depen- 
dence of plate thickness, h, on the'characteristic width, 
W, or characteristic length, L, of the particles comes 
into play. Assuming that h is sufficiently small for the 
43 impacting particle, a situation which can be checked 
and will be considered through example cases exam- 
ined later, we proceed as follows. Although the expan- 
sion waves make a measurement of U, an extremely dif- 
ficult affair, in the scheme proposed here, u„ can be 
50 measured in an easier manner. When S, reaches the 
plate rear surface, an expansion wave (FIG. 8 ) is re- 
flected back into the compressed plate material. 

Now, three separate cases must be considered and 
these are: 

Case I — Ppo > Pio 

Case II ?po Pto 

Case III ~ ppo < Pto 

Depending on which of these cases is present, different 
60 .processes occur in the plate during impact. 

Case I Pt/a ^ Pto _ . 

For this case, let us start by developing an inequality 
involving the plate compressed material velocity, 

65 and the impacting particle velocity Vp. For > p,o 
and equation ( 6 ) we have, examining FIG. 6 


Measurement of U, 
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u,i> Up, 


( 12 ) 


Therefore, substituting inequality (12) into equation 
(7), we have the following inequality. 


which may be rearranged to the following forms 
2u„ > Vp or u„ > Vppi 


( 14 ) 

Again, FIG. 8 illustrates the situation present in the 
plate after the shockwave, S,, arrival at the plate rear 
surface. Because of the zero pressure boundary condi- 
tion at the plate surface, an expansion wave is propa- 
gated back into the compressed plate material decom- 
pressing and accelerating this material giving it an addi- 
tional velocity component u , t - As the plate material al- 
ready has velocity u,,, the maximum velocity of mate- 
rial emanating from the plate rear surface (i.e., the 
maximum velocity of the spray, W/,) is given as follows: 

«/. = u„ + U,r 

( 15 ) 

It is well known that for solids, u,r can be taken as 


U,r = u„ 


( 16 ) 

with only a small error present for moderate compres- 
sions. Therefore, using this result the substitution of 
(16) into (15) and rearranging for u,, yields 

2u„ = Uf, or u„ = Uf„2. 

( 17 ) 

Now using equation (17) and inequality ( 14) an impor- 
tant observation may be made concerning the relative 
magnitudes of the maximum spray velocity, and the 
impacting particle velocity as shown in the following 
expression: 

Uf, = 2 m „ > Vp 

( 18 ) 

which yields the inequality 

Ufy > Vp. 

( 19 ) 

Inequality (19) says that the maximum spray velocity 
for the target material is greater than the impacting 
particle velocity. In other words, the plate material 
after impact disengages from the particle material. If 
the time, T,, it takes this plate material to traverse a 
known void region of width W is measured then Uf, 
would be calculated from 

Uf.= WIT, 



Once Uf, is measured u« can be obtained through equa-*^ 
tion (17). To further our understanding of (^impact 
situation discussed here for ppp > p,o conside^pte fol- 
lowing graphical interpretation of the impact with ref- 
5 erence to FIG. 9. 

The intersection of the plate and particle Hugoni^s 
gives the pressure P„ = Pp, and compressed material ve- 
locity u„ generated by the impact. A decompression of 
the material by the expansion wave, FIG. 8, from the 
10 plate rear surface is along an isentrope. If the isentrope 
is approximated by the material Hugoniot, which is 
what is implied by equation ( 1 6), then simply reflecting 
the plate material Hugoniot through the plate-particle 
Hugoniot’s intersection will yield the free surface vc- 
15 locity, i.e., the maximum spray velocity Uf, achieved 
through a shock compression and subsequent isentro- 
pic expansion. That Uf, is twice u„ is obvious from this 
procedure. FIG. 10 is a schematic of a device which al- 
lows the measurement of Uf, to be made, and FIG. 1 1 
20 is a time trace of the light level which might impinge on 
the photodiode of FIG. 10. 

When the particle penetrates the first finite thickness 
plate, FIG. 10 (a), the photodiode can record the light 
emission Fp, shown on FIG. 1 1. When the maximum ve- 
25 locity plate material having velocity «/,> Vp impacts the 
backstop [FIG. 10(I>)), a second flash F, shown on FIG. 

1 1 would be recorded. Division of the width of the void 
region W by T„ the time interval between Fp, and F,, 
yields Uf, by equation (20). From this, u„ can be calcu- 
30 lated by equation (17) and the density of the particle 
determined. The subsequent arrival time of the particle 
at the backstop would provide a second flash, Fp, at the 
time Tp later than the first flash, Fp,. Division of Tp into 
W should yield for a particle more massive than the 
33 plate the particle velocity through 


Vp=-W!Tp. 

This second velocity measurement provides an impor- 
tant method of checking the operation of the device be- 
cause Vp is known. 

What has been discussed until now is the ideal opera- 
tion of the device when a particle sufficiently large im- 
pacts a finite thickness plate. As there is no way to dis- 
criminate against particles smaller than some mass par- 
tide impacting the device, and producing spurious 
reading, the device must measure or take data which 
will allow detection of smaller than necessary particle 
impacts so that we may reject or accept with caution, 
for bounding purposes only, such data. Several exam- 
ples will now be considered of nonideal impacts, that 
is, particles less than some specific necessary size for 
ideal operation of the device. We shall now see how the 
second velocity reading (i.e., measurement of Tp and 
knowing Vp) can be used to verify the correctness of 
^ the first measurement, T,. Suppose readings of T, and 
Tp are such that in the following examples we have 


Example 1 


60 


a. Uf, < Vp 

b. Tp > Wluf, > WIVp 

Obviously this reading for the free surface velocity can 
be rejected because if ppp > p,o, Uf, > Vp in contrast to 
(a) and no use of (b) is actually required. 


65 


Example 2 


a- Uf, > Fp 

b. Tp > WIVp > Wluf,. 


( 20 ) 
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This example shows a finer use of the known impact- 
ing particle velocity, Vp, through the measurement of 
Tp. We cannot reject the data using (a) because for 
> the inequality (a) is correct. But from (b) we note 
that Tp > W/Vp, which says that the particle was 
slowed during its passage through the plate. This would 
indicate that the particle was not massive enough such 
that one portion of S, reached the plate rear surface un- 
attenuated to produce the maximum possible U/,. We, 
therefore, can reject this data or more intelligently use 
it to establish a lower bound on the particle density. 
CTase II Ppo Pto 

For this case we shall develop an equality involving 
the plate compressed material velocity, u,i, and the im- 
pacting particle velocity Vp. For p^ = p,o and equation 
(6), we have for any material from FIG. 6 

u„ = U^l 

( 21 ) 

substituting equation (21) into equation (7) we have 


10 

Example 1 

((ft ^ . 

By knowing Vp and the fact that Uf, must equal Vp for 
5 Ppo = Pto we may immediately reject this data. 

Case III Ppo < Pia 

For this case let us again develop an inequality in- 
volving the compressed material velocity and impact- 
ing particle. As before, this inequality will be used to 
to gain insight into the events which occur in the plate 
after impact. For p^ < p,o and equation (6) we have, 
examining FIG. 6, 

U,i < Upi. 

•5 ( 25 ) 

Substituting (12) into (7) we have the following in- 
equality 

20 “« ~ 

( 26 ) 


u,i = Vp- u,i 


which may be rearranged to the following forms, 
< 22 ) 25 2uu < Vp or u„ < Vpl2. 


which after rearrangement yields 


( 27 ) 


2u,i = Vp or u,i = Fp/2. 


Now we have already noted that the maximum spray 
velocity which can be reached through a single shock- 
30 wave and expansion wave is 


From equations (17) and (23) we may make the fol- Uf, = 2u,i 

lowing connection between the maximum spray veloc- 
ity and the impacting particle velocity. 

Uf, = 2uii = Vpi Uf, = Vp. 35 equations (27) and (28) indicate that 


( 24 ) 


Uf, = 2m,( < Vp or Uf, < Vp, 


( 28 ) 


Equation (23) says that the maximum spray velocity 
is equal to the impacting particle velocity. FIG. 12 
shows graphically the impact situation between a parti- 
cle and plate of the same density. 

The intersection of the plate and particle Hugoniots 
gives the pressure Pn = Ppi and the compressed material 
velocity u,i generated by the impact. A decompression 
of the material by the expansion wave, FIG. 8, from the 
plate rear surface is along an isentrope. By approximat- 
ing the plate isentrope with its Hugoniot which is the 
same as the impacting particle Hugoniot, we have at 
zero pressure that the plate material achieves a final ve- 
locity Uf„ which is the same as the impacting particle 
velocity V„. 

Using the device shown in FIG. 10, the photodiode 
would record a light emission time trace shown in FIG. 
13. When the particle penetrates the first thin plate an 
impact flash ¥p, will be sensed by the photodiode. 
When the thin plate material and particle impact the 
second plate or backstop a second impact flash F,p will 
be generated. Division of T,p from FIG. 13 into W will 
give the value of u/, which it should be remembered is 

Vp. 

Again we have considered the ideal operation of a 
device under the conditions Ppo = Pio where Uf, = Vp. 
However, if the projectile is not massive enough so that 
the thin plate achieves the impacting projectile veloc- 
ity, we will have the following nonideal example of op- 
eration for the device. 


( 29 ) 

Therefore equation (29) indicates that a thin plate of 
thickness h, and Pw > Ppo has the capability of slowing 
an infinitely large particle. Obviously, this situation is 
not possible, and what must happen is that the thin 
plate material assumes the velocity of a more massive 
particle. The process by which the plate material as- 
sumes the final terminal velocity Vp is by successive 
shock waves and expansion waves. 

Consider FIG. 14 and a massive particle impacting a 
thin plate. Immediately after impact the state in the 
^ compressed particle and plate is given by the intersec- 
tion of the particle and plate Hugoniots at P,( and u,|. 
With the arrival of the shockwave S< at the plate rear 
surface, a reflected expansion wave occurs. If the plate 
Hugoniot is used to approximate the plate isentrope, 
^ the plate material relaxes to zero pressure and a veloc- 
' ity 2u,i. The particle material would relax back to zero 
pressure and velocity Vp. However, from equation ( 27 ) 
we know that 2un < Vp. 

What we have now is a plate having velocity 2u,i 
being impacted again by a particle having velocity Vp. 
On FIG. 14 such an impact may be represented by 
starting the plate Hugoniot at pressure zero, and veloc- 
ity equal to 2 m„. The intersection of this Hugoniot with 
the particle Hugoniot gives the new pressure, P«, and 
material velocity 2 m„ -I- induced in the plate by a 
second shockwave. Reflection of this shockwave as an 
expansion wave from the plate rear siirface will generr 
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ate a second reflected expansion wave. Again approxi- 
mating the plate isentrope by the plate Hugoniot, that 
is reflecting the plate Hugoniot through P« and 2u„ -(- 
we see that the expansion wave again releases the 
material pressure back to P=0 but the material veloc- 
ity is now 2(u„ -1- «ra). The particle at this time also re- 
laxes to P = O but velocity V, where 

Pp > 2(u„ + Ua). 

(MO 

In the same manner as has been done before, we may 
consider a third impact between a particle moving at 
velocity Vp and a plate having velocity of 2(u„ + u„). 
From such successive shocks and expansion waves, we 
see that the thin plate asymptotically approaches the 
particle velocity. 

For Ppo < pt„ then Uf, = Fp is the proper conclusion 
for the maximum spray velocity for a particle more 
massive than the amount of plate material removed. 

During the operation of the device shown in FIG. 10, 
the light level impinging at the photodiode should be 
similar to that depicted in FIG. 15. 

The first flash Fp, would be the impact flash seen at 
the rear surface of the first thin sheet. The second flash 
F,p would be the impact flash due to plate and particle 
material impacting the backstop. Now, if the particle is 
sufficiently massive as compared to the amount of ma- 
terial removed from the plate division of T,p, the time 
between Fp, and F,p into W should yield 

. Fp = WIT,, 

( 31 ) 

As previously done, we have discussed the ideal opera- 
tion of the device for a projectile sufficiently massive 
that the system operated as planned. To consider how 
the system will allow us to reject data produced by the 
smaller than required particle, we need only consider 
one example. 

Example 1 

u„ = WIT,, <V, 

We see we can immediately reject the data because for 
ppo < Pio only if Uf, = V, can we be sure p„ s p„ as Case 
I, example 1, and Case II, example 1, will also produce 
a Uf, < V, if the particle is less in magnitude than the 
required particle even though p„ a p„. 

For Case I impacts discussed, the system could be 
used to determine the particle density. However, for 
Case III impacts, the same system (remembering Case 
II impacts) will only indicate that p„ s p„ but not the 
value of p,o- Certainly this situation indicates that one 
would select plate materials of the lowest possible den- 
sity for use in the device. Then, certainly, even knowing 
Vpo s Pio would be information having engineering and 
scientific significance. 

In Case III discussed heretofore, the desirability of 
using the lowest density material available to construct 
the first thin plate, for a metal plate either beryllium or 
magnesium, was noted. To illustrate, assume beryllium 
is used as the thin plate material. What will be of con- 
cern is the maximum beryllium material spray velocity 
when particles having greater or equal density impact. 
In addition, for impacting particles having densities less 
than the beryllium plate, the beryllium plate velocity 
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after two shock and expansion wave transits will be de- 
termined. 

Using the graphical techniques developed to answer 
these concerns, the beryllium Hugoniot is plotted start- 
5 ing at pressure and material velocity equal to zero on 
FIG. 16. Selecting an arbitrary impact velocity, say V, 
= 0.7 cm/microsec. to have a definite example, the 
Hugoniots for gold ( 19.24 gm/cm®), iron (7.8 gm/cm^), 
iodine (4.93 gm/cm’), aluminum (2.79 gm/cm''), beryl- 
10 lium (1.85 gm/cm^), magnesium (1.725 gm/cm '), car- 
bon dioxide (1.54 gm/cm'), and water (1.0 gm/cm') 
are plotted in a reflected manner starting at pressure 
equal to zero and material velocity equal to V„. Either 
by examination of the intersection of the various mate- 
1 5 rial Hugoniots with the beryllium Hugoniot or by calcu- 
lation using equation (9), the following values of u,„ 
Ua, and therefore, Uf„ were obtained and are displayed 
in the Table of FIG. 17. 

When measurements are received from the device 
20 such that Uf, > y„ the values in the table may be used 
to determine or bound the impacting particle density. 
For example, if Uf, = 0.98 cm/microsec. then the parti- 
cle must be of a density which is that of iron. For Uf, = 
0.78 cm/microsec., the iodine and aluminum values 
25 form an upper and lower bound for Uf, and, therefore, 
the particles density is less than that of iodine but 
greater than that of aluminum. To obtain a better value 
of the particle’s density additional Hugoniots for mate- 
rials having densities greater than aluminum but less 
39 than iodine may be used. 

Considering the second item of interest, it can be 
noted from FIG. 17 that, for materials down to 1 
gm/cm®, in two-shock and expansion wave transits 
Uf,iy,> 0.081. 

33 Although the invention has been described relative to 
specific embodiments thereof, it is not so limited and 
the described specifics are for illustration of the princi- 
pal and are not intended to serve as limitations thereof. 
Thus, modifications and variations in the specific em- 
bodiments described will be readily apparent to those 
skilled in the art in the light of the above teachings. It 
is therefore to be understood that within the scope of 
the appended claims, the invention may be practiced 
otherwise than as specifically described herein. 

3 What is claimed as new and desired to be secured by 
Letters Patent of the United States is; 

1. A method to determine the density of unknown 
material particles wherein the closing velocity is known 
between the impacting particles and a plate of known 
material comprising: 

determining the material velocity that would be pro- 
duced in the plate by the unknown material parti- 
cle of known closing velocity; 

determining the material velocity that would be pro- 
duced in the plate by various known material parti- 
cles having the same known closing velocity upon 
impact with the plate; 

deriving the unknown particle material density by 
making use of the facts that material velocity in the 
plate is ordered on impacting particle density alone 
with all other variables held constant and that simi- 
lar density impacting particles having the same 
density produce the same material velocity in the 
^5 same material plate; 

by comparing the material velocity said unknown 
particle produced in the known material plate and 
the material velocities that were produced in the 
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plate by the various known different material parti- ' 
cles upon impact; 

obtaining a coincidence between the unknown parti- 
cle material velocity in the plate and a known parti- 
cle material velocity in the plate; and, 5 

measuring the plate material by; 
measuring the velocity of material exiting from the 
plate opposite the particle impact side by recording 
(with a photomultiplier or photon tube) the initial 
particle plate impact flash; 10 

starting a time counter; 

recording the exiting plate material impact flash on 
a backstop (with a photomultiplier or photo tube), 
stopping the time counter, dividing the recorded 
time into the known distance between plate and is 
backstop to obtain the exiting material velocity and 
, dividing the exiting material velocity by two to ob- 
tain the plate material velocity. 

2. A method to determine the density of unknown 
material particles wherein the closing velocity is known 20 
between the impacting particles and a plate of known 
material comprising: 

determining the shockwave velocity that would be 
produced in the plate by the unknown material par- 
ticle of known closing velocity; 25 

determining the shockwave velocity that would be 
produced in the plate by various known material 
particles having the same known closing velocity 
upon impact with the plate; 
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deriving the unknown particle material density by 
comparing the shockwave velocity said unknown 
particle produced in the known material plate and 
-the shockwave velocities that were produced in the 
plate by the various known different material parti- 
cles upon impact; 

obtaining a coincidence of the unknown particle' 
shockwave velocity in the plate with a known parti- 
cle shockwave velocity in the plate; 
making use of the facts that shockwave velocity in 
the plate is ordered on impacting particle density 
alone with all other variables held constant and 
that similar density impacting particles having the 
same density produce the same shockwave velocity 
in the same plate; and, 
measuring the shockwave velocity by: 
starting a time counter upon recording the impact 
flash; 

recording the impact flash on the plate front surface; 

determining the time the shockwave reaches the 
plate rear surface by recording the impact flash of 
plate and argon gas; and 

stopping the time counter, to thereby derive the 
shockwave velocity by the time it takes to travel 
from the front to rear surface of a known finite 
thickness plate. 

* ^: * * * 
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